Thermodynamic and spectral properties of compressed Ce calculated by the merger 
of the local density approximation and dynamical mean field theory 
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We have calculated thermodynamic and spectral properties of Ce metal over a wide range of 
volume and temperature, including the effects of 4/ electron correlations, by the merger of the 
local density approximation and dynamical mean field theory (DMFT). The DMFT equations are 
solved using the quantum Monte Carlo technique supplemented by the more approximate Hubbard 
I and Hartree Fock methods. At large volume we find Hubbard split spectra, the associated local 
moment, and an entropy consistent with degeneracy in the moment direction. On compression 
through the volume range of the observed 7-0: transition, an Abrikosov-Suhl resonance begins to 
grow rapidly in the 4/ spectra at the Fermi level, a corresponding peak develops in the specific heat, 
and the entropy drops rapidly in the presence of a persistent, although somewhat reduced local 
moment. Our parameter-free spectra agree well with experiment at the a- and 7-Ce volumes, and a 
region of negative curvature in the correlation energy leads to a shallowness in the low-temperature 
total energy over this volume range which is consistent with the 7-a transition. As measured 
by the double occupancy, we find a noticeable decrease in correlation on compression across the 
transition; however, even at the smallest volumes considered, Ce remains strongly correlated with 
residual Hubbard bands to either side of a dominant Fermi-level structure. These characteristics 
are discussed in light of current theories for the volume collapse transition in Ce. 
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PACS numbers: 71.27.+a, 71.20.Eh, 75.20.Hr 



I. INTRODUCTION 



A number of rare earth metals undergo pressure in- 
duced first order phase transitions with unusually large 
volume changes of 9-15% (for reviews see Refs. |l|,||). 
Of these transitions the isostructural 7-0; transition in Ce 
has received the most attention. □ ft was discovered first, 
has the largest volume change (15% at room tempera- 
ture) , and may also be accessed entirely at ambient pres- 
sure (or in vacuum) by changing the temperature, thus, 
for example, allowing thorough spectroscopic investiga- 
tion of both phases. The results of such photoemission 
and Bremsstrahlung studiesO show a dramatic transfer of 
spectral weight to the Fermi energy and the development 
of a large peak with its center of gravity slightly above the 
Fermi energy when going from the 7- to the a-Ce phase. 
Similarly, the optical conductivity is higher in the a phase 
where the frequency dependent scattering rate is charac- 
teristic for a Femii liquid behavior with an effective mass 
of about 20 me .El Also the magnetic susceptibility and its 
temperature dependence change from a Curie- Weiss like 
behavior in the phase to a Pauli paramagnetic behavior 
in the a phase.Q Despite these dramatic differences, the 
number of 4/ electrons does not change significantly and 
is close to one across the 7-0; phase transition line which 
ends in a critical point at T = 600 ± 50Kjj above which 
the 7- and a-Ce phases become indistinguishable. 

Notwithstanding the considerable attention, there re- 
mains continued disagreement about the nature of the 
transition and the a phase. In general, it is believed 
that the transition is driven by changes in the 4/ elec- 



tron correlations, though some alternative theories have 
been proposed. Two recent examples of the latter as- 
sume some kind of symmetry breaking in the a phase: 
Eliashberg and CapellmannEl argue that a-Ce has a sym- 
metry broken distorted structure mainly based on the 
observation thai, the a phase shows large changes of the 
compressibility ;E| Nikolaev and MichelO propose (hidden) 
quadrupolar ordering. In these theories a critical end- 
point is impossible because of the symmetry breaking, 
and there must be, at least, a second order phase transi- 
tion line above 600 K, which disamees with the com- 
mon interpretation of the experiments 

The first theory of electronic origin tp-describe the 7-a 
transition was the promotional modeM which assumed 
a change in the electronic configuration from 4f^(spd)^ 
to Af^{spd)'^. However, it was soon ruled out by ex- 
periment which did not reveal any major change in the 
number of 4/ electrons. Also band structure calculations 
found about one 4/ electron per Ce atom in the a phase, 
leading Johansson to propose a Mott transition (MT) 
scenario.llil Similar to the MT of the Hubbard model, 
the 4/ electrons are considered to be localized in the 7 
phase and to be itinerant in the a phase, with this reduc- 
tion in the degree of 4/-electron correlation being caused 
by the decrease in the ratio of Coulomb interaction to 
kinetic energy under pressure. In a subsequent analy- 
sis based on these ideas, Johansson et altB employed a 
standard local density approximation (LDA) calculation 
for the spdf electrons in the a phase, while treating the 
4/ electrons as localized 4/^ moments decoupled from 
LDA spd bands in the 7 phase. Evidence for the MT 
scenario to be correct is taken from the considerable sue- 
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cess of LDA calculations and their generalized gradient 
improvements for the structural and vol un i ie .dependence 
of the total energy of a-Ce-like phases.EJtJ Addiiioiial 
support appears to come fror norbi t ally polarizcdIljIlj'tZI 
and self-interaction correctedOHEj LDA modifications 
which obtain transitions in Ce and Pr at about the right 
pressures. Also LDA+U calculations have been reported 
for one or both Ce phasesEJ^El 

This MT model appears in conflict with the Kondo vol- 
ume collapse (|K.VC) scenario of Allen and Martincj and 
Lavagna et alB3 which is based on the Anderson impu- 
rity model. Both pictures agree that, at the experimental 
temperatures, the larger volume 7 phase is strongly cor- 
related (localized), has Hubbard split 4/ spectra, and 
exhibits a Hund's rule 4/^ moment as reflected in the 
observed Curie- Weiss magnetic susceptibility. But, while 
the MT scenario then envisages a rather abrupt transi- 
tion on compression to a weakly correlated (itinerant) 
a phase, in which the Hubbard split bands have coa- 
lesced together near the Fermi level and the 4/^ moments 
are lost, the KVC picture assumes continued strong cor- 
relation in the a phase with Kondo screening of the 
4/ moments by the valence electrons. The signature 
of this Kondo screening is a peak in the 4/ spectra at 
the Fermi level, the Abrikosov-Suhl resonance, which 
lies between the remaining Hubbardpsplit spectral weight 
characteristic of the local moment J3 Assuming a rapid 
volume dependence of the Kondo temperature, there is 
Kondo screening and an Abrikosov-Suhl resonance in the 
a phase but not in the 7 phase since the Kondo temper- 
ature of the a phase is above and that of the 7 phase be- 
low the typical experimental temperatures. The strong 
volume dependence of the Kondo temperature leads to 
a region of negative curvature in the free energy and, 
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similar to a 



thence, a first order transition as diet, 
mon (Maxwell) tangent construction '^'^'^^ 
vapor-liquid transition. While the KVC model provides 
a genuine many-body calculation for Ce, it nevertheless 
incorporates only two bands, which has prompted pre- 
vious jattempts at introducing more orbital realism from 
LDA.BEa 

On the basis of the underlying Hubbard and periodic 
Anderson models, it has been argued recently that the 
MT and KVC scenarios are rather similar when the im- 
portant many-body effects are taken into account .£3 That 
is, the behavior of the local moment at the MT of the 
Hubbard model is not so abrupt, nor is the appearance 
of a three peak structure in the density of states unique to 
the periodic Anderson model. This fact can be obscured 
by the use of static mean-field approximations (including 
LDA and its modifications) especially when describing 
the a phase. 

A new approach to describe Ce including both or- 
bital realism and electronic comjlaticn effects is now 
available with the recent mergeiOHc^ 
namical mean field theory (DMFT),E 
has been employed by Zolfl et alS3 who used the non- 
crossing approximation (NCA) to solve the DMFT equa- 
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tions in order to calculate the spectra, Kondo temper- 
atures, and susceptibilities for a- and 7-Ce. Indepen- 
dently, we treated the DMFT equations with the more 
rigorous Quantum Monte CarlotS (QMC) simulations 
and reported, as first results of the present effort, evi- 
dence for a Ce volume collapse in the total LDA+DMFT 
energy which—coincides with dramatic changes in the 
4/ spectrumO A similar transition was also |described 
earlier in LDA+DMFT calculations for PuBj In all 
three cases, the / spectra showed Abrikosov-Suhl reso- 
nances lying in between residual Hubbard splitting for 
the smaller- volume, less-correlated a phases, in contrast 
to the LDA results mentioned above which only obtain 
the Fermi-level structure. Related behavior is also ob- 
served for the Mott transitioruip V2O3, which was stud- 
ied recently by LDA+DMFT. H 

In the present work we extend Ref. ^ to lower teruper-. 
atures, complement it with Hubbard-I calculations,E3LZl 
and calculate the volume-dependence of additional phys- 
ical quantities including the entropy, specific heat, total 
spectrum, orbital occupation, and the magnetic moment. 
In Section II, the LDA+DMFT technique is briefly de- 
scribed along with the Hubbard-I approximation and a 
new and faster implementation of the QMC treatment 
which is subsequently validated against established ap- 
proaches. In Section III, thermodynamic results, i.e., the 
energy, specific heat, entropy, and free energy, are pre- 
sented over a wide range of volume and temperature and 
the signatures for the a-7 transition are discussed. We 
present the volume- and temperature-dependence of the 
4/- and the valence spd-spectrum and compare to exper- 
iment in Section IV. The 4/ occupation, local magnetic 
moment, and related quantities are given in Section V. 
Finally, the results of this paper are summarized and dis- 
cussed in Section VI. 



II. THEORETICAL METHODS 

The results in this paper have been obtained by the 
LDA+DMFT method, that is by the merger of the lo- 
cal density approximation (LDA) and dynamical mean- 
field theory (DlMFT) which was recently introduced by 
Anisimov et alBJl and Lichtenstein and Katsnelsonc3 (for 
an introduction see Ref. ^9|). The starting point of this 
method is a conventional LDA band structure calcula- 
tion. Since electronic correlations are only treated at a 
mean field level within LDA, the most important term 
for electronic correlations, i.e., the local Coulomb inter- 
action, is added explicitly. This defines a multi-band 
many-body problem which is solved by DMFT. To solve 
the DMFT equations, we employ two different implemen- 
tations of the quantum,Monte Carlo (QMC) technique as 
well as the Hubbard-lE3CZl (H-I) approximation. This sec- 
tion describes the relevant computational details of our 
calculations. 
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A. LDA+DMFT approach 

Scalar-relatpastic, linear muffin-tin orbital LDA 
calculationsLJO were performed for face centered cubic, 
(fcc) Ce over a grid of volumes as described elsewhere.cl 
The associated (6s, 6p, 5d, if) one-electron Hamiltonians 
define 16 x 16 matrices H^^j^, after shifting the 4/ site en- 
ergies to avoid double counting the Coulomb interaction 
Uf between 4/ electrons. The latter is explicitly taken 
into account in the full second-quantized Hamiltonian for 
the electrons, 

H = ^ (H^DaO^)) Im.l'm' c\ilma^kl' m' a 

U..lrnA'm'.(T 



(1) 



i,mcr,m'(T' 



Here, k are Brillouin zone vectors, i are lattice sites, 
Im denote the angular momentum, a is the spin quan- 



tum number, Ui/^ 



, and the prime signi- 



fies ma ^ m'a' . The many-body Hamiltonian Eq. (|l|) 
has no free parameters since we employed constrained- 
occupation LDA calculations to determine Uf and the 
4/ site energy shift for all volumes considered (see Fig. 5 
of Ref. Hlfor the values). We did not take into account the 
spin-orbit interaction which has a rather small impact on 
LDA results for Ce, and also neglected the intra-atomic 
exchange interaction which has only an effect if there are 
more than one 4/-electrons on a Ce atom. 

The DMFT maps the lattice problem Eq. (|l|) onto the 
self-consistent solution of the Dyson equation 

Gkitiu) =[tiul + ^il - i?V(k) - S(zc^)//] , (2) 

and a seven-orbital (auxiliary) impurity problem defined 
by the bath Green function 



number electrons described by Eq. (|^) is n/-|-n„ = 4 per 
Ce site. Here, the number of 4/ electrons n/, and simi- 
larly the number of valence (i.e., spd) electrons, may 
be obtained from the lattice Green function 



(4) 



where T is the temperature and — {2n+l)T:T are the 
Matsubara frequencies. To obtain the physically relevant 
Green function G{uj)^ we employ the maximum entropy 
methodcj for the analytic continuation to real frequencies 

LO. 

In principle the LDA and DMFT parts of the calcu- 
lation should be mutually self-consistent, with DMFT 
changes in orbital occupations (especially Uf) feeding 
back into a ncw-J/Lp^ (k) and Uf , as argued by Savrasov 
and coworkers .E3 Certainly the constrained occupation 
calculations used to fix Uf and the 4/ site energy in 
-^LDA(k) should not be impacted, as they are intended 
to be valid over the range 0<nf<2. These calculations 
provide what are in effect the screened Coulomb ener- 
gies for 0, 1, and 2 / electrons per site, which covers this 
range according to what fraction of|i;he sites are at one or 
another of the various occupations. □ However, differences 
between the DMFT n f and the LDA n / could, if the for- 
mer were fed back into the LDA, change the position of 
4/ level slightly, and with that the extension of the 4/ 
wavefunction, and thus the /-valence hybridization. It is 
simply not known at this point if such effects are impor- 
tant, although we note that DMFT(QMC) and LDA-like 
(see Sec. [ID) solutions of Eq. (|l|) generally yield values 
of n/ within 10% of one another. The additional cost on 
top of the already very expensive LDA-I-DMFT(QMC) 
method also makes such additional self-consistency im- 
practical in the present case. 



g{ioj)-'^\^Y.^r{G^{iu^)If}] +E(*c.). (3) 



Here / is the unit matrix. If = [5if5iif5mm'] projects 
onto the seven /-orbitals, ji is the chemical potential, 
Tr denotes the trace over the orbital matrix, and N is 
the number of k points [N = 2048 for T < 0.4 eV and 
A;' = 256 for T > 0.4 eV). Within the LDA, there is a mi- 
nor crystal-field splitting of the seven 4/ orbitals. How- 
ever, in Eq. (^) we average over the seven 4/ orbitals, i.e., 
we treat them as degenerate in the auxiliary impurity 
problem. Consequently the DMFT self-energy is diago- 
nal Yi{iuj)If, at least in the paramagnetic phase studied. 
The impurity problem is solved with one of the meth- 
ods described in the following two sections and generates 
a self-energy T,(iLu). This self-energy gives a new Green 
function in Eq. (H) and thus a new impurity problem and 
so on, iterating to self-consistency (for more details see 
Refs. ^ and 29). In this self-consistency cycle, the chem- 
ical potential /i of Eq. (0) is adjusted so that the total 



B. Hubbard-I approximation 

In the large-volume limit where intersite hybridization 
vanishes, the auxiliary impurity problem is simply the 
isolated atom, i.e., ^(iw) = l/{iuj+^ — ef) where e/ is the 
4/ site energy. In this limit the exact self-energy is known 
and may, at finite-^fokimes, be used as the Hubbard-I (H- 
I) approximationcaO: 



Wj(Mat,T) 



3=^ 



-KA^at-(j-l)C//' 



14T^G^*(zw„)e*'^"°^, 



(5) 
(6) 

(7) 



where £/ has been absorbed into /iat, which is set at 
each iteration in such a way that n^* of Eq. (|^) equals 
the current iif of Eq. (^. The positive weights Wj for 
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transitions between j — I and j electrons are given by 



= [jVj + il5-j)v,_,]/iuY,^i)^ (8) 



1=0 



where Vj are Boltzmann weights for having j electrons 
on the atom 



14! 



' j!(14-j)! 



exp [_{ij(j-l)C//-jMat}/r].(9) 



Our DMFT(H-I) procedure is in fact also correct at all 
volumes in the high-temperature limit. Noting that the 
Wj^s sum to unity, one can see that 

S-*M = UfY,jv,/Y,v, = -ufUf, (10) 



3=0 3=0 



since we always choose /J,at so that n^* —n-f. This is the 
paramagnetic Hartree-Fock value, which is also the cor- 
rect high-temperature limit since the w„ oc T, and only 
the high-frequency tail of the self-energy is of importance. 



C. QMC simulations 

Our main approach to solve the DMFT impurity prob- 
lem is the numerical QMC technique. We use two imple- 
mentations which differ mainly by the Fourier transfor- 
mation between the Matsubara frequency representation 
employed in the Dyson equation Eq. (|^) and the imagi- 
nary time representation employed for the QMC simula- 
tion of the impurity problem Eq. (||). Within QMC, the 
imaginary time interval [0,(3] (/3 — ^/T) is discretized 
into A Trotter slices of size Ar = /3/A. Since there are 
91 Ising fields per time slice, the number of time slices 
which are computationally manageable in the QMC is se- 
riously restricted. Thus, if one employs a discrete Fourier 
transformation between Qiiujn) at a finite number of A 
Matsubara frequencies and G{ti), ti = IAt, . . . , AAr, the 
resulting Green function oscillates considerably around 
the correct G{A To overcome this shortcoming, Ulmke 
and coworkeraiil suggested using a smoothing procedure 
which replaces 5(iw„) — s- t/(iw„), after calculating the 
auxiliary Giiujn) via Eq. (|^), where 



1 - e-xp[-AT/g{iujn)] 



(11) 



It is g{iLUn) that is Fourier transformed to imaginary 
time g{Ti), and once the QMC simulations of the An- 
derson impurity problem have yielded the output G(t;), 
the process is reversed: The Fourier transform of G(t;), 
G(iujn), yields G{iLUn) from the inverse of Eq. |ll| The 
new self-energy is then I](iw„) — g{iujn)^^ — G{iujn)^^. 

This approach generates smooth Green functions G{ti) 
and reproduces the correct At — » limit. We use it in 
one implementation of the QMC algorithm, referred to 



as QMCi in the following. Other approaches employed 
in the literature are to fit splines to G(t;) and, thus, 
to use more suoport points than A to do the Fourier 
transformationEll or to extend the Matsubara frequency 
sums by employing theJlerated perturbation theory re- 
sult at high frequencies. l3 Most results of our paper were 
obtained by yet another QMC implementation (QMC2) 
which uses a different way to Fourier transform and which 
seems to be less sensitive to statistical noise. As this 
modification is new, it is described in some detail in Sec- 
tion 



II C 1 and validated in Section [I C 2 



1. Modified QMC implementation 
In the implementation QMC2, we use a constrained fit 
G(r)=5]u;J,(T). (12) 



to the output QMC impurity Green function G{ti), in 
order to accomplish the Fourier transform to G(ia;„) for 
n = -^N^, . . . , ^N^ - 1 with > A. The basis func- 
tions arc /i(r) = —e^^^'^/{e~^'^' + l) and have Fourier 
transforms fi(iuj) — 1/ [iuj — Si). At real frequencies, 
Eq. (12) corresponds to a set of (5- functions with differ- 
ent spectral weights wi, and is capable in the limit of an 
infinite set of basis functions of reproducing any given 
spectrum. In contrast to a spline-fit where every fit- 
coefficient is determined by the local behavior in an imag- 
inary time interval, in our approach every fit-coefficient 
is determined by the local behavior in frequency space. 

The constraints to the fit Eq. (O) are Wi > 0, G(0''") 
is precisely the QMC value, G(0^-F G(/3"^ ) = -1, and 
^G(0+) + ^G(/3-) =g2, where p- = /3-0+ and g™ is 
the (iw)^'" high-frequency moment of G{iuj). For the last 
constraint, g2 is obtained from the relation G^^{iuj) = 
Q~^[itjj) — Y^iiijj) which implies g2 = 52 + Sq, where these 
are the indicated moments of G{ioj), Q[iuj), and 'E{iuj), 
respectively. Here, 52 is known as G is input to the QMC, 
and we take sq = I](zCi; = oo) = (13/14)n/[/y with tLf, — 
14[l-|-G(r = 0+)] for the present paramagnetic case.E3 

Typically we use grids of A/4 equally spaced Si, and 
optimize the agreement with the QMC data as a function 
of the centroid and width of these grids, in each case sys- 
tematically eliminating basis functions for a given grid 
which would otherwise yield negative Wi. Because the 
QMC expense increases as A'^, we are forced to execute 
fewer Monte Carlo sweeps for the largest A's, and the 
statistics become less good than for smaller A's. How- 
ever, the constraint Wi>0 still seems to provide a sensi- 
ble interpolation through the statistical noise, although 
this has the consequence that the number of surviving 
positive Wi increases more slowly than A. Nonetheless, 
we see a systematic evolution as a function of A and 
extrapolations A—^oo agree with large- volume and high- 
temperature limits (Hubbard-I) and the QMCi (see Sec- 
tion |IIC2D. Note that while we find the fit Eq. (O) to 
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be very useful for functional behavior along the imagi- 
nary time and frequency axes, and for integral quantities 
such as Uf and the total energy, it is not useful in prac- 
tice for directly obtaining real frequency behavior in the 
presence of typical QMC statistical uncertainties. The 
maximum entropy method is far superior here as iti'olds 
these uncertainties into calculation of the spectra.E2l 

In order to accelerate the convergence of our 
DMFT(QMC2) method we carry out cheap iterations 
on the constant part of the self-energy in between each 
expensive QMC iteration. That|-is, we subtract a con- 
stant Hartree-Fock contributionEj from the QMC self 
energy: AE(icj) = Y.^^'^iiuj) - {13/U)nJ^'^Uf where 
14[1 + G(t = 0+)]. Following every QMC cycle. 



QMC _ 



then, one has 'E{iLu) = (13/14)n/L// -f AI](iw) in Eq. (H) 
which is iterated to self-consistency with nj from Eq. (jj), 
while keeping AS fixed. The resultant values of Uf and 
nj^'^ agree within statistical uncertainties. These un- 
certainties can be significantly smaller for nf than for 
rij'^'^ at the smallest volumes. 

We find G(t) to converge quickly as a function of 
QMC2 iteration for all r at small volume, and for r close 
to and /3 at large volume. For intermediate t at large 
volume and low temperature, however, where G(t) is 
generally quite small, convergence appears to result from 
the average of frequent small values of G(r) with occa- 
sional large values as the Ising configurations are sam- 
pled, with the large-volume atomic limit approached by 
the latter becoming statistically unimportant. 

In order to improve the statistics given this behavior, 
we have chosen to include sweeps from all previous QMC 
iterations (excluding warm up sweeps) along with the 
new sweeps in Gf™{Ti) in arriving at the QMC2 result 
for iteration i: G,(tO = [Gf'''^{Ti) + (i - 1)G,_i(tO] A- 
Note that the warm-up sweeps themselves are already 
started with a reasonable self-energy, such as a converged 
DMFT(H-I) resuh or a DMFT(QMC2) result for another 
A. We have tested this treatment at both small and large 
volumes by starting anew at i = 1 from the converged 
DMFT(QMC2) self-energy, and have found agreement 
with the previous results to within statistical uncertain- 
ties. 

We used 10, 000 sweeps per QMC iteration for A = 80, 
decreasing systematically to 1,000 for A = 256, and car- 
ried out from 20 to over 100 QMC iterations for each 
T, V point. At small V even at T = 0.054 eV we found 
the DMFT(QMC2) energy to settle down generally after 
a few QMC2 iterations to maximal excursions of about 
±0.02 eV (±0.05 eV) for A = 80 (256), with the root- 
mean-square uncertainties much smaller. Such benign 
behavior extends to increasingly large volumes at higher 
T, where these DMFT(QMC2) results begin to agree 
closely with DMFT(II-I). At low temperature, the scat- 
ter in our measurements as a function of iteration grows 
as volume is increased, especially in the transition region 
and beyond; however, the Trotter corrections also dimin- 
ish here so there is less need for larger A. 



Finally, we turn to the issue of performing the Fourier 
transform from imaginary time to Matsubara frequen- 
cies. The virtue of the fit Eq. ( |T2| ) is that it decouples 
A and allowing manageable QMC costs (smaller A) 
and yet accurate kinetic energies (larger N^). Most of our 
DMFT(QMC2) calculations took iV^^256 for T> 0.054 
eV and = 512 for T = 0.027 eV. In the course of this 
work we realized that there is a volume dependence to 
the error in the kinetic energy from the Matsubara cutoff, 
and while the Ni^ = 256 choice at T = 0.054 eV leads to a 
small 0.04 eV error in the vicinity of the transition, it be- 
comes more significant, 0.11 eV, at the smallest volumes 
considered. Since our DMFT(H-I) and DMFT(QMC2) 
codes have identical kinetic energy treatment, we used 
the former to correct the present DMFT(QMC2) results 
to effective values of four times those just noted, 
which should give better than 0.01 eV accuracy at all 
volumes. We verified this by selected DMFT(QMC2) 
tests with the larger N^. Note that this kinetic energy 
treatment includes (and the cited errors reflect) an ap- 
proximate evaluation of the full infinite Matsubara sum. 
Specifically, we approximate the high-frequency behavior 
of a quantity F(iLo) by Ft){iLjj)=wi/ (iLJ—ei)-\-W2/ (iw— £2), 
with parameters chosen to reproduce its l/(ia;)™ mo- 
ments for m = 1-4. Then we approximate the infinite 
Matsubara sum on F{iuj) by the analytic result for the 
infinite sum over Fq plus the finite sum over the difference 



2. Validation 

H ere we validate the new faster QMC2 algorithm of 
Sec. II C 1 , used for much of the low-temperature ther- 



modynamic results in this paper, against QMCi which 
employs the Ulmke-smoothing. Such validation involves 
extrapolation to the limits iV^, A — > 00, where the QMCi 
approach should provide exact results. Errors which van- 
ish in these limits include those arising from truncation 
of Matsubara sums (finite N^), and from the Trotter ap- 
proximation (finite A). 

Figure |^ compares the kinetic DMFT energy (see Sec- 
tion II D| for details of its calculation) obtained by QMC2 
and QMCi as a function of At = /9/A, at a temper- 
ature T = 0.54 eV and atomic volume V = 16.8 A^. 
(Note that the Ar dependence is largest at small vol- 
umes, as we shall discuss further.) The line with open 
circles shows the QMCi results with Matsubara sums 
truncated after = A frequencies under the applica- 
tion of Ulmke's smoothing procedure Eq. (|Tl|). Those 
with squares and open triangles show the results when 
these sums are extended to A^^^ = 00 using the Hartree- 
Fock (HF) Green function at high Matsubara frequen- 
cies; that is, using Eq. with E ^ Shf = (13/14)n/[// 
for — AttT, . . . , oo.eET In the first case (squares) the 
current chemical potential and Uf were used to define 
Shf- In the second case (open triangles) the whole pro- 
cedure was made self-consistent: From Uf, we calculate 
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FIG. 1; Extrapolation of the kinetic energy at r = 

0.54 eV and l/ = 16.8A^ using the QMC2 (fiUed circles) and 
the QMCi algorithm (open circles, squares, and triangles; 
differences are due to whether and how Hartree-Fock results 
for the high-frequency tails of self-energy are included, see 
text). The lines show the extrapolations through the QMC 
data yielding £;kin(Ar = 0) = 49.888 ±0.003 eV (filled circles) 
and 49.853 ± 0.022 eV (open triangles), both with a mixed 
quadratic and cubic fit; and 49.944 ± 0.271 eV (open circles) 
and 49.713 ± 0.305 eV (open squares), both with a linear fit. 
The results agree within twice the above standard deviation 
and, thus, validate, the QMC2 algorithm. The inset shows 
the two upper curves (filled circles and open triangles) as a 
function of Ar^ over an expanded At range. For the QMC2 
(filled circles) it also compares the mixed quadratic/cubic fit 
(solid line) with a purely quadratic fit to the data points which 
fulfiU ATUf/2 < 0.4 (dotted line). 



Y,{iuj) — {13/ 14:)n fU f + AY,{iuj) for all frequencies at a 
fixed AT,{iLo) which is defined in the previous section. 
This T,{iuj) yields a new n/ via Eq. (Q), and so on until 
convergence. As can be seen, the dependence on At is 
greatly reduced, as is a lso the case for the QMC2 im- 
plementation of Section II C 1 (filled circles) which also 
uses this self-consistent treatment of the HF part of the 
the self-energy. To avoid a large At error, the large fre- 
quency part of the self-energy and especially the con- 
stant HF part is important to the energy, and must be 
self-consistently correct. 

The inset of Figure |^ shows the top two curves in an ex- 
panded scale versus At^, which indicates 0.035 eV agree- 
ment between the two QMC methods in the limit At^O. 
We believe that effects of Matsubara cutoff are largely 
eliminated here and that Trotter errors predominate in 
these curves, which are expected to be of leading order 
At^ , at least for the QMC calculation of lattice |Bi©dels 
without the DMFT self-consistency complicationH'Ca In 
order to keep the Trotter errors under control, it is cecom- 
mended that At be constrained to At[/j/2 < 1,EJ and 
we have done so in this work. In fact at small volumes 
where the Trotter corrections are the largest, we would 



find the need to use three terms over this full range, with 
both a + bAr^ + cAt^ and a + bAr + cAt^ providing rea- 
sonable fits to the energy. We find the ratio c/b to be sig- 
nificantly smaller for the first choice, which is consistent 
with the expectations of a leading At^ dependence. In 
our DMFT(QMC2) calculations we have therefore cho- 
sen to use the two term fit a -I- 6At^, however, over a 
reduced range. The dotted line in the inset of Fig. |l|, for 
example, suggests that At/7 //2 < 0.4 might be a reason- 
able range for this fit, given Uf = 5.05 eV for the volume 
in the figure. The two-term fit also makes more sense 
in the volume range of the transition, where the Trotter 
corrections are smaller, but there is also more scatter as 
a function of At. 



In DMFT(QMC2) calculations for the whole volume 
grid we have used at least 1,2, and 3 At values for tem- 
peratures greater than, equal to, and less then 0.544 eV, 
respectively. In the first case it is easy to take At so 
small that really no extrapolation is needed, or maybe 
one other value as a spot check at the smallest volume. 
At T = 0.054eV (632 K) we used At[///2 = 0.417, 0.334, 
and 0.209.EI Our calculations at T = 0.027 eV (316 K) 
were limited by expense to systematically larger values, 
At[///2 = 0.667, 0.477, and 0.334, so that extrapolations 
to At = are more uncertain. Even the smallest At here, 
which corresponds to A = 320, leads to a At^ value that 
is 2.6 times larger than its counterpart at T = 0.054 eV. 
Fortunately, we see every indication that our electronic 
Hamiltonian is already very close to its low-temperature 
limit by T = 0.054 eV (632 K), as these total energies 
agree with those at r = 0.027 eV within their error bars 
at the same finite At values. The At— >0 extrapolations 
are more benign for Uf and d which also agree well for 
the two temperatures. Accordingly, we do not display the 
r = 0.027eV results in this paper, but do comment on 
the agreement between the two temperatures as specific 
quantities are presented. 



We have alluded earlier to the fact that the Trotter 
approximation errors get larger at smaller volume in the 
present work. This makes sense as these are related to 
the commutator of the kinetic and potential energies, 
and should thus depend on the size of the hybridization, 
which gets larger as volume is reduced. We find no dis- 
cernible dependence of the energy on At for volumes in 
the 7 phase of Ce for the range of A investigated, but 
that for the smaller volume a phase, we find dE/dAr^ 
to become significant and to increase in magnitude with 
decreasing V. Since the 7-0; transition is intrinsically 
related to the growing importance of hybridization ver- 
sus the Coulomb interaction as volume is decreased, this 
behavior is perhaps not surprising, although the effect 
turns on rather abruptly, appearing as almost another 
signature of the transition. Similar beii^vior has been 
seen for the periodic Anderson model.! 
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D. Calculation of the LDA+DMFT energy 

There are several possible expressions for the DMFT 
total energy per site, depending on whether the potential 
energy is obtained using the self energy 

i^DMFT = § E [{^fLDA(k) + ^^{l^n)} 

X Gk(*c^n)]e^""°\ (13) 
or from a thermal average of the interaction in Eq. (|l|) 

^DMFT = ^ E [^LDA(k)Gk(lC^„)] 6^""°^ +Ufd. 

(14) 

In the latter expression, 

I ' 

E {f^ifmaflifm'a') (15) 



2N 



i ma^m'a' 



is a generalization of the one-band double occupation for 
multi-band models, which may be calculated directly in 
the QMC presuming the site average is given by the as- 
sociated impurity problem. If we were using the exact, k- 
dependent self-energy in these equations, Eq. (13) would, 
be equivalent to Eq. (|l^ and to the Galitskii-MigdalE^ 
expression for the total energy. We find Eq. (14) to be far 
superior in the present LDA-I-DMFT(QMC) calculations 
in regard to low-temperature stability and agreement 
with known limits, possibly not surprising in that it takes 
a thermal expectation of the true Coulomb interaction for 
the problem. We use Eq. (Ill) for the LDA+DMFT(H- 



I) energy in preference to Eq. (14) with a purely atomic 
calculation of d. However, it should be noted that for 
vanishing intersite hybridization at large volume, as well 
as at high temperatures, the H-I self-energy is exact, and 
indeed we find agreement between QMC and H-I results 
for -Edmft in these limits. 

To evaluate the total LDA-I-DMFT energy Etot {T) in- 
cluding all core and outer electrons, we add a correction 
to the paramagnetic all-electron LDA energy i?LDA(r) 

^tot(T) = ^LDa(T) + ^DMFt(T) - S,nLDA(r) , (16) 

which consists of the DMFT energy i?DMFT(7') less an 
LDA-like solution of the same many-body model Hamil- 
tonian Eq. (|l]), thus "model LDA" or £;mLDA(T). The 
latter is achieved by a self-consistent solution of Eqs. @) 
and (^ for n/ taking a self-energy SmLDA = Uf{nf — 
From this, the kinetic energy is calculated by the first 
term of Eq. (|lj) and the potential energy by ^Ufnf{nf~ 
1). Note that while all of thase expressions are explic- 
itly temperature dependent ,Cj the present calculations 
are electronic only and do not attempt to add lattice- 
vibrational contributions. Estimates of these contribu- 
tions are similar for the a- and 7-Ce phases, howevec, 
and appear to have little impact on the phase diagram.E3 



Finally, one virtue of the Hamiltonians Eq. (P is that 
it is possible to reach high-temperature limits where the 
entropy is precisely known. One may then calculate the 
entropy from the DMFT energy 



'DMFT 



^y,^ rfi?DMFT(r) ^ (17) 



dT' 



where Soo = kB[M\nM - n\nn - (M-n) In (M-n)] = 
12.057 fee, for M = 32 states and n — A electrons per site 
for Ce. 



III. THERMODYNAMICS 

In this section we consider thermodynamic properties 
of Ce, more specifically the energy, specific heat, entropy, 
and free energy, over a wide range of volume V and tem- 
perature T. Intercomparison of the Hartree-Fock (HF), 
DMFT(H-I), and DMFT(QMC) methods to solve the ef- 
fective LDA Hamiltonian Eq. (m here serves to validate 
all three calculations in limits where they should and do 
give the same answers, and also to point out shortcomings 
of the more approximate techniques elsewhere. Then we 
turn specifically to the a-7 transition in Ce, and use the 
total energy, Eq. (16), and entropy, Eq. ([l7|), to present 
evidence for the volume collapse transition. 



A. Global behavior 

Figure || shows the correlation energy of the effective 
LDA Hamiltonian Eq. (|^), defined as the energy E of 
Eq. (|l]) less the paramagnetic HF result iJpMHF for the 
same Hamiltonian. Results for polarized HF, DMFT(H- 
n^^and DMFT(QMC) as obtained from Eqs. and 
( p^ ) in the last two cases, respectively, are compared in 
this manner for an extended range of atomic volumes 
at five temperatures. The polarized HF solutions as- 
sume ferromagnetic spin order, and display both spin 
and orbital polarization, with one band depressed be- 
low the Fermi level and the other thirteen lying above. 
These HF calculations (dash-dot curves) are seen to give 
good energies at large volume and low-temperature in 
comparison to the DMFT(QMC) (solid lines with data 
points), as has beem observed previously for the Ander- 
son HamiltonianJlZlo Thus, the polarized Hartree-Fock 
solution and other polarized stati£_Hiean-field methods 
such as orbitalUtjpjilarized LDA,y£aO self-interaction- 
corrected LDAjtZlllal^ and LDA-|-UEj£]I can be expected 
to give good low-temperature total energies in the strong 
coupling limit. As the atomic volume is reduced, the dif- 
ference between polarized and paramagnetic HF energy, 
-BpoiHF — £'PMHF, becomes positive near 22 A'', and the 
HF solution has a transition from the polarized to the 
paramagnetic solution, where all fourteen bands have co- 
alesced together above but slightly overlapping the Fermi 
level. The highest two temperatures in Fig. ^ lie above 
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FIG. 2: Correlation energy, i.e., the difference between the 
LDA+DMFT and the paramagnetic HF (PMHF) energy, as a 
function of volume at five temperatures. At large volumes, the 
LDA+DMFT(QMC) energy agrees with the polarized HF and 
the Hubbard-I (H-I) solutions. But the LDA+DMFT(QMC) 
energy breaks away from the polarized HF energy for decreas- 
ing volume, leading to a region of negative curvature in the 
vicinity of the experimentally observed Q-7 transition (indi- 
cated by arrows) at low temperature. 



the critical point for this transition, and so there is no 
polarized HF solution. 

Turning to the Hubbard-I approximation, which be- 
comes exact in the atomic limit, it is no surprise that the 
DMFT(H-I) results (short-dashed curves) should agree 
well with the DMFT(QMC) energies at large volume for 
all temperatures. This approximation is also exact in the 
high temperature limit, as may be seen from Fig. |^, where 
there is also increasingly good agreement at high temper- 
ature for all volumes considered here. The agreement be- 
tween the two distinct DMFT calculations in these limits 
provides a test of the reliability of both approaches used 
here. 

A direct view of the temperature dependence is given in 
Fig. ^ where the energy E of Eq. (|l|) is plotted versus T 
for an atomic volume V of 46 A^. At this relatively large 
volume, the DMFT(QMC) and DMFT(H-I) resuhs agree 
closely and smoothly interpolate between the polarized 
HF energy at low temperatures and the paramagnetic 
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FIG. 3: Energy (upper figure) and entropy (lower figure) 
of the LDA-I-DMFT Hamiltonian Eq. (|l|) vs. temperature at 
y = 46 A^. At this relatively large volume, the DMFT(QMC) 
and DMFT(H-I) energies agree with each other and, at lower 
temperatures, also with the polarized Hartree-Fock solution. 
However, the entropy of the latter is completely wrong since 
the 14-fold degeneracy of the local magnetic moment is disre- 
garded. 



HF result at high temperature (above about 15 eV, not 
shown). There is no temperature-induced transition in 
the DMFT results here, in contrast to the unphysical 
transition from the paramagnetic to the polarized phase 
within HF at T'--^ 1 eV. This transition is a shortcoming of 
the paramagnetic HF phase in which double-occupations 
of /-electrons on the same Ce site can not be avoided 
such that the paramagnetic (interaction) energy is too 
high. 

Additional insight is provided by the corresponding en- 
tropy in Fig. The DMFT(H-l) entropy approaches 
/cb ln(14) at low temperature, which is effectively the de- 
generacy of the Hund's rules magnetic moment /cb ln(2 J-t- 
1), where without intra-atomic exchange and spin orbit 
interaction we get the full 14-fold degeneracy of the / 
level rather than the proper 6-fold degeneracy for J = 5/2. 
At still lower temperatures, crystal field effects are known 
to reduce the entropy.CJ 

Figure ^ illustrates two important aspects in which HF 
and more rigorous techniques differ. First, the HF transi- 
tion at about 1 eV corresponds to simultaneous moment 
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formation and magnetic ordering. In contrast, the two 
processes are distinct in more rigorous treatments, with 
the moment formation occurring in a continuous fashion 
at higher temperatures, culminating in the low-T plateau 
in Fig. I^b, with the onset of magnetic order (if it occurs) 
coming at yet lower temperatures off the scale of the plot. 
Second, polarized HF gives good low-T energies at large 
volumes because one of the Hund's rules multiplet states 
will be a single Slater determinant. However, its broken 
symmetry mistreats the entropy at lower temperatures, 
giving zero instead of, e.g., /cBln(2J+l) for nf — 1 in 
the atomic limit, so that the finite-T thermodynamics 
are incorrect. 



2.5 



B. 



Transition 



We now consider thermodynamic evidence for the a- 
7 transition in Ce. While the QMC error bars restrict 
us from making a quantitative prediction, we argue that 
the present results do suggest the transition. Evidence 
is already apparent in Fig. ^, where the DMFT(QMC) 
correlation energy is seen to bend away from the polar- 
ized HF result as temperature is lowered, leading to a 
region of negative curvature in the vicinity of the ob- 
served transition (arrows). As the other terms (-Blda 
and -EpMHF ~ ^-iiiLDa) contributing to the total energy 
Eq. (|l^) all have positive curvature throughout the range 
considered in this work, this correlation contribution is 
then the only candidate to create a region of negative 
bulk modulus in the low-temperature total energy, i.e., a 
thermodynamic instability, and thence a first order phase 
transition given by the Maxwell common tangent. 

Figure ^ shows total energies Eq. (|6|) for the 
DMFT(QMC) and polarized HF methods at the three 
lowest temperatures of Fig. ^ The region of negative cur- 
vature just noted in the correlation energy is seen to cause 
a substantial depression of the DMFT(QMC) total ener- 
gies (solid curves with symbols) away from the polarized 
HF results (dashed curves) below 35 A'^, which is most 
pronounced at the lowest temperature, r = 0.054 eV. The 
consequent shallowness in the DMFT(QMC) curve at 
this temperature persists over the observed range of the 
two-phase region (arrows), although statistical uncertain- 
ties preclude any claim of seeing negative curvature. The 
slope is also consistent with a —0.6 GPa pressure (long- 
dashed line),_which is the extrapolated transition pres- 
sure at r = 0£3 We suggest in fact that these T = 0.054 eV 
(632 K) total electronic energies are already close to the 
low-T limit. Both the DMFT(H-I) and HF energies at 
this temperature differ by less than 0.006 eV from cor- 
responding results at half this temperature, throughout 
the volume range in Fig. |. Our DMFT(QMC) calcu- 
lations at T = 0.027 eV (316 K) are also consistent with 
this conclusion, as discussed in Sec. HC2| , 

That the electronic contribution to the total energy 
might be close to its low temperature limit below about 
600 K is also consistent with the analysis of the a-j tran- 
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FIG. 4: Total LDA-HDMFT(QMC) and polarized HF en- 
ergy as a function of volume at three temperatures. While 
the polarized HF energy has one pronounced minimum in 
the 7-Ce phase, the negative curvature of the correlation en- 
ergy of Fig. ^ results in the development of a side structure 
(r = 0.14eV), and finally a shallowness (r = 0.054eV), which 
is consistent with the observed a-7 transition (arrows) within 
our error bars. These results are also consistent with the ex- 
perimental pressure given by the negative slope of the dashed 
line. 



sition by Johansson et a/.,@ who attribute the temper- 
ature dependence of the transition pressure primarily to 
the difference in entropy, which is zero and /cb In (2 J-f 1) 
for the a and 7 phases, respectively. That is, for tem- 
peratures larger than hqik the Kondo temperature and 
the crystal-field splittingca in the 7 phase, yet still fairly 
low (say 200-600 K), the temperature dependence of the 
7-phase free energy may be dominated by the linear term 
—fee In (2 J-|- 1)T arising from a plateau such as in Fig. ^d, 
while presumably the total energies (both a and 7) are 
closer to the low-T limit due to their faster T^ depen- 
dence. 

We have calculated both the DMFT(QMC) specific 
heat C{V,T) and entropy S{V,T) for the effective Ce 
LDA Hamiltonian Eq. (§). We first calculated DMFT(H- 
I) energies Eq. ( |l3| ) on a logarithmic temperature grid 
up to the high-T limit (^ 10"^ eV) where the entropy is 
known to be 12.057 fee- As noted earlier, the DMFT(H- 
I) method is correct at high temperatures, and indeed 
the DMFT(QMC) energies obtained via Eq. (|l|) closely 
approach the H-I results as T is increased, e.g., lying 
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above by only 0.024 and then 0.004 eV at T = 5.4 
and 13.6 eV, respectively, for V = 16.8 A^. We there- 
fore fit the difference between the QMC and H-I ener- 
gies at eight temperatures from 0.054 to 5.4 eV to the 
form a + ^^^n/i^ + nA/T'^), n = 1-3, which has a 
behavior at low temperatures, and is benign at high 
temperatures. These smoothed and interpolated differ- 
ences were added to the DMFT(H-I) energies to cre- 
ate a fine grid of "DMFT(QMC)" energies from which 
C(V,T) = dE{V,T)/dT\v was calculated by numerical 
differentiation, and S{V^ T) by integration down from the 
high-T limit according to Eq. (p7|). Note that while the 
finite nature of Eq. (|l|) is unphysical at very high tem- 
peratures, these results are nonetheless entirely meaning- 
ful at more modest temperatures where the omitted core 
and higher-lying valence states will be frozen out, e.g., 
below ^ 3 eV near the a-7 transition, given a spectrum 
of Eq. that extends to nearly 30 eV above the Fermi 
level in that volume range. 

The challenging need for accurate energy derivatives, 
as well as the sensitivity of Eq. (^7|) to the lowest temper- 
atures given the 1/T factor, requires a stringent conver- 
gence criterion for the kinetic-energy Matsubara sums. 
Otherwise we observe unphysical negative low-T limits 
of the entropy for V < 25 A'^. We have also constrained 
the fits to smooth out the value of this low-T limit as a 
function of volume over this same range. In all cases it 
is to be emphasized that the fits give excellent represen- 
tation of features in -Bdmft(qmc) (^) - -£^dmft(h-i)(7'), 
ranging in size from 0.1 to 0.24 eV upon decreasing the 
volume from V ^ 35 to 25 A^, and are well within the 
±0.03 eV error bars in the data. The same fits were used 
to obtain both C{V,T) and S{V,T). 

The temperature dependence of the DMFT(QMC) 
specific heat is shown in Fig. ^ at six volumes. The 
most significant feature is the appearance of the low tem- 
perature peak in the range T = 0.1-0.2 eV, which co- 
incides precisely with growth of the quasiparticle peak 
or Abriksov-Suhl resonance at the Fermi level in the 4/ 
spectra, as will be seen in the next section. Analogous 
behavipc has been discussed for the one-band Hubbard 
model.Eil The low temperature peak in the specific heat is 
just barely discernible at the 7-phase volume of 34 A'^ in 
Fig. ^, has become rather prominent by 29 A^, which is 
slightly larger than the a-phase volume, and then contin- 
ues to broaden and shift to higher temperatures as vol- 
ume is further reduced. The broad peak near 1 eV which 
appears at all volumes is due both to the 4/ charge fluc- 
tuations, and also to spd- valence to 4/ excitations, given 
that n f increases by ~ 20% on raising the temperature to 
1.4 eV. Note also in regard to the charge fluctuations that 
the peak in C(T) should occur at significantly smaller T 
than the Coulomb repulsion [// ~ 6 eV, as may be seen 
in the case of the half-filled one-band Hubbard model. 
Here, the specific heat peak occurs at r = 0.208 U in the 
absence of hopping t = 0, and the lOjCation of the peak is 
also depressed by the band widthH 

The volume dependence of our DMFT(QMC) entropy 




Temperature T (eV) 



FIG. 5: Specific heat as a function of temperature for dif- 
ferent volumes (oflF-set as indicated). At smaller volumes, an 
additional low-energy peak develops, coinciding with the for- 
mation of an Abrikosov-Suhl resonance (see Fig. below). 



is shown in Fig. ^ for six temperatures. The rapid 
increase in the entropy over the a-7 transition (28.2- 
34.4 A^) is due precisely to the low temperature peak 
peak in C(r), which contributes to the entropy via its 
weighted area / dTC{T)/T. Thus, at large volumes 
where the 4/ spectral weight is Hubbard-split with no 
contribution at the Fermi level, the low-T entropy is 
pinned at fcBln(2J-|-l) (ignoring effects of crystal field 
at yet lower T). Then, as the volume is reduced, the 
quasiparticle peak begins to grow at the Fermi level, the 
weighted area of its associated heat capacity peak reduces 
the low-T entropy below A:b ln(2J-|-l) via Eq. (|l^. The 
physical interpretation is of course that the degeneracy 
associated with the 2J-I-1 directions of the Hund's rules 
moment disappears as this moment is either screened or 
collapses on reducing the volume. 

For completeness, we conclude this section by provid- 
ing the free energy F = Etot - ST in Fig. @, although 
the uncertain errors in the entropy, and the fact that the 
large-y, low-T value is 50% too large [taking into ac- 
count the spin-orbit coupling will give ln(6) instead of 
/cBln(14)]. Given that the electronic total energy Etot 
is near its low-T limit by T = 0.054 eV, we consider 
that curve as "T = 0", and then include it again as 
F = Etot- ST for T = 0.054 eV. The error bars on aU 
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FIG. 6: Entropy as a function of volumes for different tem- 
peratures. In the vicinity of the a-'y transition (28.2-34.3 A''), 
the entropy increases rapidly. 



curves are just from the energy. The slopes of the two 
straight lines give the experimental transition pressures 
at T = and 0.054 eV, and arrows mark the observed 
boundaries of the a-7 transition at room temperature. 
The essential conclusion of Fig. |^ is that these results 
are consistent with experiment, though stronger claims 
are precluded by the statistical uncertainties. Nonethe- 
less, the results of this section which we find compelling 
are the way in which EtotiV) systematically develops a 
shallowness in the vicinity of the a-7 transition as tem- 
perature is lowered, and the structure in the specific heat 
and entropy. 



IV. SPECTRA 

In this section, we discuss the spectral changes through 
the a-7 transition. To obtain the physical spectrum 
A{ijj) — — ilmG((jj), one has to analytically continue 
the QMC data from the imaginary time (Matsubara fre- 
quency) representation to real frequencies oj: 



G{t) 



duj ^7 tA(uj) 

00 -1- L 



(18) 



> 




V(A') 

FIG. 7: Free energy as a function of volume at three tem- 
peratures, compared to lines whose negative slopes give the 
experimental a-7 transition pressures at T = (solid line) and 
0.054 eV (dashed line). Given the statistical uncertainties, the 
results are consistent with experiment and show that a shift 
of the Qf-7 transition volumes is primarily due to the entropy. 



regime. To deal with this ill-conditioned problem which is 
particularly cumbersome in the presence of the statistical 
QMC error, we employ the maximum entropy method.Cj 
When interpreting the results later on, we have to keep in 
mind, however, that there is a significant error at larger 
frequencies which tends to smear out fine features such 
that, e.g., inner structures of Hubbard bands are not nec- 



essarily resolved. In section IV A , we present the spectra 
of the f- and valence-electrons of fee Ce as a function 
of volume and discuss the changes at the a-7 transition. 
The spectra obtained are compared to photocmission and 



Bremsstrahlung experiments in Section IV B 



A. Change of the spectra at the a-7 transition 



As one can see in Eq. (|18[), the values of A{u;) at large 
(positive or negative) frequencies affect G{t) only weakly 
because the integral kernel is exponentially small in this 



In Section III we noted a region of negative curvature in 
the correlation energy at volumes consistent with the ex- 
perimental a- and 7-volumes, leading to a shallowness in 
the total energy and suggesting a first order phase tran- 
sition at lower temperatures. To further elucidate the 
nature of the ongoing changes, we study the evolution of 
the /-electron spectrum as a function of volume for fee 
Ce at r = 0.054 eV (632 K) in Fig. |. This temperature 
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is close to the critical endpoint (T = 600 ± 50 K) at whicb. 
the first order a-7 transition disappears experimentally.u 
From the continuous evolution of the energy versus vol- 
ume curves, we expect, however, similar changes above 
the critical endpoint, which are not yet strong enough to 
cause a first order phase transition. At a very small vol- 
ume, V = 20A^, most of the spectral weight is seen to be 
in a big quasi particle peak or Abrikosov-Suhl resonance 
at the Fermi energy, but some spectral weight has already 
been transferred to side structures which would be inter- 
preted as upper and lower Hubbard bands in a Hubbard 
model. Moving closer to the a-7 transition (between 28.2 
and 34.4 at room temperature), the a-Ce-like spectra 
at y = 29 shows this three peak structure to become 
more pronounced with a sharp Abrikosov-Suhl resonance 
and well-separated Hubbard bands. The spectral weight 
of the Abrikosov-Suhl resonance is further reduced and 
smeared out when going to the 7-phase {V = 34 A'^) and 
finally disappears at large volumes {V = 46 A^), at least 
at T = 0.054 eV. 

Altogether, we observe, as a function of volume, the 
crossover from a structure which differs only slightly from 
a one-peak structure, to a three peak structure, and fi- 
nally to a two peak structure. The physical interpre- 
tation is that the /-electrons are somewhat correlated 
at low volumes, where the large quasiparticle peak above 
the Fermi energy resembles (to a first approximation) the 
one-peak structure of the uncorrelated one-particle the- 
ory or the LDA. At larger volumes, the system is highly 
correlated, there is a magnetic moment imposed by the 
electrons in the lower Hubbard band, but the /-electrons 
at the Fermi energy are still itinerant. Finally at the 
largest volumes, the /-electrons are localized and the lo- 
cal magnetic moment is fully developed. Here, the most 
dramatic change of the weight of the quasiparticle peak 
coincides with the observed region of negative curvature 
in the correlation energy. We thus conclude that the 
drastic reduction of the weight of the quasiparticle peak 
is related to the energetic changes in the correlation and 
total energies which are consistent with the first order 
a-7 transition. 

These features and also the three peak (Kondo-likc) 
structure of the a and 7r^hases agree with the Kondo 
volume collapse scenario.E3 On the other hand, many- 
body calculations show that the behavior of the Anderson 
and Hubbard models — paradiems for the Kondo volume 
coUapsecj and Mott transitionO scenarios, respectively 
— are remarkably similar in regard toJieir spectra and 
other properties at finite temperatures Ej One important 
difference, however, is the absence of spectral weight at 
the Fermi level in the "large volume" phase ofptJpe Mott- 
Hubbard transition, as for example in V2 03j^ in con- 
trast to the reduced but still extant jSpectral weight in 
our 7-Ce results and the experiment Jj which is a more 
Kondo-like feature. 

In Fig. 1^, we compare the 4/-spectrum in the a and 
7 phases to results at higher temperatures (T = 0.14eV) 
from Ref. Isl. Most notably, the Abrikosov-Suhl reso- 
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FIG. 8: Evolution of the 4/ electron spectrum with volume 
a.t T — 632 K; off-set as indicated. When going from small to 
large volume, the weight of the central Abrikosov-Suhl reso- 
nance decreases and practically fades away at the Q-7 transi- 
tion from = 29 to 34A^ . The residual weight around the 
Fermi energy at = 34A'^ indicates a smeared out Abrikosov- 
Suhl resonance as is to be expected if the Kondo temperature 
of 7-Ce is below T = 632 K. 



nance in the a-phase {V = 29 A'^) becomes much sharper 
when going from r = 0.14eV to 0.054 eV. The reason for 
this is that the Abrikosov-Suhl resonance is smeared out 
thermally at T = 0.14eV (1580 K) since this temperature 
is comparable to the Kondo temperature, which we esti- 
mate to be 0.18 eV (2100 K) from the full width at half 
maximum 
K, see Ref. 



LDA+DMFT(NCA) calculations yield 1000 
This Kondo temperature is only a crude 



32 



estimate which might also be somewhat reduced if the 
spin orbit coupling, which splits off states from the Fermi 
energy, is taken into account. Nonetheless, it reasonably 
agrees with experimental estimates of Tk = 945 K and 
1800-200Q.K for the Kondo temperature obtained_from 
electronicp and high-energy neutron spectroscopyLd, re- 
spectively. In contrast, the peak in the 7 phase remains 
smeared out such that one would assume a Kondo tem- 
perature lower than 0.054 eV (632 K); the experimental 
estimates are Tk = 95K (Ref. |) and 60 K (Ref. |l]). 
The changes in the rest of the spectrum are much less 
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FIG. 9: 4/ electron spectrum for a-Ce {V = 29A^) and 7- 
Ce {V = 34A^) at two temperatures (T = 632 K: solid line; 
T = 1580 K: dashed line). The Abrikosov-Suhl resonance of 
a-Ce is smeared out when increasing the temperature from 
T — 632 to 1580 K, indicating that the Kondo temperature 
is in between. 



FIG. 10: Evolution of the spd electron spectrum with volume 
at T = 632 K; off-set as indicated. Note the wider energy 
window in comparison to Figs. ^ and ^. The main effect to 
be seen is the decrease of the bandwidth upon increasing the 
volume. 



B. Comparison to experiment 



dramatic. The position of the upper Hubbard band is 
fixed while the lower Hubbard band, which has a very 
small spectral weight, moves closer to the Fermi energy 
upon decreasing the temperature. 

While the /-electrons undergo a transition from itin- 
erant character at low volumes with a quasiparticle reso- 
nance at the Fermi energy, to localized character at larger 
volumes without such a resonance, the (spd) valence elec- 
trons remain metallic at all volumes. This can be seen 
in Fig. ^ which shows the valence spectral function 
A{oj) averaged over the spc?-orbitals. It is finite at the 
Fermi energy for all volumes, such that Ce is always a 
metal. The biggest change in the spectrum is the de- 
creasing valence bandwidth when increasing the volume, 
which is simply due to the reduced overlap of the valence 
orbitals as inter-atomic distances increase. The effect 
of electronic correlations is less obvious. But, one can 
note a dip in the valence spectrum in the vicinity of the 
Fermi energy which is to be expected to coincide with the 
Abrikosov-Suhl resonance in the /-spectrum. This dip is 
most pronounced at lower volumes where the /-electron 
Abrikosov-Suhl resonance has most spectral weight. 



The LDA-HDMFT(QMC) calculation of fee Ce sug- 
gests a volume collapse approximately at the experimen- 
tal volumes. To further test whether this theory actu- 
ally describes fee Ce, we now compare our a- an|d-|7-Ce 
spectra with photoemission spectroscopy (PES)E3_,and 
Bremsstrahlung isochromatic spectroscopy (BIS)H To 
this e nd, we combined the / and spd spectra of Section 
IV A with areas normalized to 14 and 18, respectively, to 
yield the full spdf density of states, and smoothed it with 
the experimental resolution of approximately 0.4 eV. 

The comparison is shown in Fig. |ll| for a and 7 
Ce. Although there are no, free parameters in our 
LDA-hDMFT(QMC) results^ the agreement between 
theory and experiment is very good. Particularly good is 
the agreement of the spectrum around the Fermi energy 
for both a and 7 Ce; this part of the spectrum consists of 
the Abrikosov-Suhl resonance of the /-electron spectrum 
and the valence spectrum. Also the position of the up- 
per and lower Hubbard bands and the relative weight of 
these peaks and the Abrikosov-Suhl resonance is correctly 
predicted by the theory. Less good is the agreement with 
respect to the width of the upper Hubbard band which is 
too narrow in our theory; the experimental upper Hub- 
bard bands extend to energies 1-2 eV higher than our 
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FIG. 11; Comparison erf the LDA+DMFT(QMC) spectra 
with experiment (circles)u. Although there are no free pa- 
rameters in the calculated spectrum, the agreement is very 
good, in particular at the Fermi energy (oj = 0). The addi- 
tional structure in the upper Hubbard band which is seen in 
the experiment is likely due to the exchange interaction which 
was neglected in our calculation. 



theory. As has been argued in Ref. this can be under- 
stood by the Hund's rules exchange couphng which has 
not been taken into account in our calculation. We justi- 
fied this by noting that the exchange couphng is only ef- 
fective if there are more than two electrons on one Ce site 
which happens only rarely. However, the excited states of 
the upper Hubbard-band correspond to just such double 
occupied states. For these, the Hund's rules coupling be- 
comes important and will split the upper Hubbard band 
into multiplets. With this shortcoming resolved, the com- 
parison to the experimental spectrum suggests that our 
LDA-hDMFT(QMC) calculation describes a and 7 Ce 
very well. 

The a and 7 spectra of mevious LDA+DMFT(NCA) 
calculations by Zolfl et alX3 are considerably different 
from ours and the experimental spectra, in particular the 
weight of the upper Hubbard bands was much higher in 
Ref. ^ The temperature of Ref. ^ is very close to ours 
(T = 580 K) and also the 4/-clcctron Coulomb interaction 
value C// is comparable, at least for the 7 phase; Zolfl. et 
al. employed a fixed value of J7/ = 6 eV whereas the con- 
strained LDA values in our calculations are C// = 5.72eV 
and 5.98 eV for a and 7 Ce, respectively. In view of this 
we tend to explain the differences, at least for 7-Ce, by 



the different method employed to solve the DMFT equa- 
tions, in particular, since the non-crossing approximation 
(NCA) is a resolvent perturbation theory for strong cou- 
pling. 



V. LOCAL MOMENT AND 4/ OCCUPATIONS 

Important additional information about the a-7 tran- 
sition and the effects of electron correlation in Ce are 
contained in the number of 4/ electrons per site n/, the 
double occupation d, and quantities derived from these 
such as the fraction of sites with ri = 0,1,2 / elec- 

trons, and the local magnetic moment. These parameters 
can discriminate between the various models, as for ex- 
ample the promotional modelll3 assumes a considerable 
change in the number of 4/ electrons at the a-^|^;ran- 
sition, in contract to the Kondo Volume coUapseEa and 
Mott transition^ scenarios which do not. The latter two 
on the other hand distinguish themselves by assuming a 
small and large change of the magnetic moment, respec- 
tively. 

Figure ^ siyes nj as a function of volume at four 
temperatures .E3 The lowest curve at T = 0.054 eV (632 
K) is already very close to the low-T limit, as our results 
at half this temperature are the same to within gener- 
ally 0.004, or at most 0.01 electrons per site. There are 
two main tendencies: With decreasing V, Uf increases 
due to the upward motion of the 6s, p levels relative to 
the 4/ level under compression; it also increases with T 
due to the thermal occupation of the large 4/ density 
of states lying above the Fermi level. Superimposed on 
this behavior is, at low temperatures, an abrupt reduc- 
tion of n/ in the observed two-phase region (marked) as 
volume is reduced, an anomaly which is annealed away 
by T = 0.5 eV similar to the case of the total energy. This 
effect leads to a number of 4/ ^ctrons close to one, rul- 
ing out the promotional modelE3 and suggesting Kondo 
physics given the sharp quasiparticle peak seen in the 
previous section. Quantitatively, we get a 4% reduction 
in nf across the two-phase region from 1.035 ± 0.017 to 
0.993 ±0.010 at r< 0.054 eV. Similar behavionis seen in 
the 10% drop from 1.014 to 0.908 of Zolfi et alH in their 
LDA-|-DMFT(NCA)calculations, and the 11% reduction 
from 0.971 ± 0.Q06 to 0.861 ± 0.015 electrons/site of Liu 
and coworkers,El who fitted a single impurity Anderson 
model to the experimental 4/ spectrum. The reason for 
the drop in is a systematic increase in the double oc- 
cupation d under compression. Since d is the potential 
energy divided by C//, the energy cost associated with 
increasing d can be ameliorated by reducing n f . 

Since there is little chance of more than doubly- 
occupied sites in Ce at low temperature, nf and d provide 
sufficient information to obtain the fractions of sites with 
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FIG. 12; Number of 4/ electrons vs. volume at four 
temperatures. At low temperatures and in the vicinity of the 
a-7 transition, n/ is very close to one. 



various integral /" occupations. 

W{f) = l-Hf+d 

w{f^) = uf -2d 



(19) 



Fig. |T3| shows our DMFT(QMC) results for these weights 
at T = 0.054 eV, which are also close to the low- 
temperature limit. At large volume one sees that each 
site nearly always has one / electron, and that empty 
or doubly occupied sites are rare, as would be expected 
for 71/ ~ 1 in the absence of significant hybridization to 
move these electrons to either f oi v {spd valence) states 
on neighboring sites. For the / electrons to begin to 
move around from one site to another in any indepen- 
dent fashion under the influence of larger // hybridiza- 
tion, or for there to be virtual charge fluctuations of the 
form f^v^ — > f'^v'^ and f^v^ — > f'^v^ due to increased 
fv hybridization, it is clear in each case that both empty 
and doubly occupied sites must become more common 
at the expense of singly occupied sites if the volume is 
reduced, as evident in Fig. |l^. Note that these changes 
are especially dramatic over the experimental two phase 
region (marked). 



The filled symbols in Fig. 13 show the impurity Ander- 
son model results of Liu et alB at the observed a- and 
7-Ce volumes; the large open symbols, the DMFT(NCA) 
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FIG. 13: Fraction of empty (/"), singly (/^), and doubly oc- 
cupied sites (/^) vs. volume as calculated by LDA-I-DMFT 
(QMC) (open symbols with lines) at T = 0.054 eV ip-pom- 
parison to LDA-I-DMFT(NCA) (large open sy mhols)B and 
impurity Anderson model results (filled symbols) While the 
DMFT results agree very well for the 7 phase, there are sig- 
nificant differences in the a phase as discussed in the text. 



results of Zolfl and coworkers.^ Our DMFT(QMC) val- 
ues are wif ),w{p),w{p) = 0.013 ± 0.019 (0.118 ± 
0.025), 0.939 ± 0.028 (0.771 ± 0.033), and 0.048 ± 0.009 
(0.111 ± 0.008) for the 7 (a) volumes, respectively. The 
two DMFT calculations agree well within these uncer- 
tainties for all three populations u'(/") at the larger 
7-phase volume, and also with the impurity Anderson 
model value for w{f^); although for the two small pop- 
ulations, they obtain w(/'') < w{£^) in reverse order to 
the values of Liu and coworkers. S The most significant 
difference at the a-Ce volume is the rather larger double 
occupancy, d = w{f'^) = 0.111 ± 0.008, obtained by our 
DMFT(QMC) calculations in comparison to smaller val- 
ues 0.044 and 0.026 obtained by the the DMFT(NCA) 
and impurity Anderson model, respectively. Temper- 
ature is unlikely to be a factor here, as we obtain nf 
and d unchanged within our error bars at half the tem- 
perature of the DMFT(QMC) results in Fig. |l|, e.g., 
d = w{p) = 0.108 ± 0.008 at T = 0.027eV (316 K). 

There are some differences between the three calcula- 
tions, however, which might account for differing wlP) 
predictions: (i) Our calculated Coulomb interaction for 
a-Ce, Uf = 5.7 eV, is slightly smaller than the Uf = 6 eV 
employed in Refs. |3^ and ^. (ii) Liu et al. employ an 
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impurity Anderson model whereas both we and Zolfl et 
al. extract a periodic Anderson type of model from the 
LDA, including /-/ hybridization. While we also deal 
with an Anderson impurity model in the course of our 
DMFT solution, this impurity model is only an auxiliary 
construction with a complicated and strongly renormal- 
ized (non-constant) hybridization, (iii) Finalbi. in con- 
trast to the DMFT(QMC), both DMFT(NCA)B and the 
1 /N approachE3 of Ref . H are based on perturbation ex- 
pansions in the hybridization strength, a quantity which 
gets larger with reduced volume. Thus, while these two 
approximations are controlled by the smallness of the hy- 
bridization strength and also by 1/N (we have 7V= 14 / 
orbitals), there are nonetheless larger corrections when 
the hybridization is increased, i.e., when going to the 
more itinerant a-Ce. Note in this context that the ra- 
tio of Coulomb interaction to an effective bandwidth de- 
termined by the total /-/ and /-valence hybridisation 
changes from 3.8 to 2.5 across the 7-0; transition.Q 

It is possible to quantify the degree of /-electron cor- 
relation by noting certain limiting values of d. A natural 
minimum is provided by the strongly correlated ground 
state of Eq. (|l|) in the atomic limit, where d is a piece- 
wise linear function of rt/, with d^dmin — max(0,nj^ — 1) 
for 71/ < 2. Similarly, d^ax = (13/28)n2 from Eq. (|l|) 
in the uncorrelated limit (riin2) = (ni)(n2), which is ap- 
proached for volume V —tQ leading to a vanishing ratio 
of Coulomb interaction to bandwidth. Fig. ^ shows a 
plot of the ratio (dmax — rf)/(dmax — rfmin) for the present 
Ce calculations at T = 0.054 eV, which reflects strong and 
weak correlation limits at 1 and 0, respectively. Note the 
polarized to paramagnetic HE transition at 20 A'^ 
for decreasing volume, and the fact that the paramag- 
netic HE result is completely uncorrelated [d k, c?max) as 
expected. The fact that the d ratio in this case is not pre- 
cisely zero is due to a small amount of orbital polarization 
arising from the fact that 4/ bands of different symme- 
try overlap the Fermi level to slightly different extent, 
whereas dmax was defined for all spin-orbital occupations 
to be 7iy/14. 

The combination of increasing d and decreasing ny 
causes a sharp decrease in correlation (delocalization) of 
the DMET(QMC) result for decreasing volume through 
the observed 7-0; transition (marked) ,, in a greement with 
tenets of the Mott-transition modeLEirEO The value of 
the DMET(QMC) d ratio is 0.76 ± 0.08 at the a vol- 
ume, combining all of the uncertainties in both d and 
n f . While this_value is certainly less correlated than the 
DMET(NeA)E3 (large open circles) and impurity Ander- 
son modefi (filled circles) predictions at 0.89 and 0.92, re- 
spectively, it is far from the kind of uncorrelated behavior 
seen in the paramagnetic HF of Fig. |lj or, presumably 
also, in the LDA. Even at the smallest volumes consid- 
ered, the DMFT(QMC) d ratio still suggests the pres- 
ence of significant correlation, which is entirely consis- 
tent with the remnant Hubbard side bands in this range 
as discussed in the previous section. Most notable in the 
DMET(H-I) curve is a glitch at about V=17A^ which is 
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FIG. 14: Double occupation ratio (dmax — rf)/(dmax — dmin) 
and local magnetic moment {m%) (triangles) as a func- 
tion of volume at T = 0.054 eV. In the former case, 
we compare the LDA+DMFT(QMC) results with our HF 
and LDA+DMFT(H-I) results as well as with the LDA-f- 
DMFT(NCA) by Zolfl et alH and the Anderson model cal- 
culations by Liu et al.U The double occupancy increases when 
going from 7- to a-Ce (experimental volumes as indicated), 
i.e., the electrons become more itinerant or less correlated. 
This effect is most pronounced in our LDA-I-DMFT(QMC) 
results; however, the d ratio is still far from the uncorrelated 
value d — dmin, i.e., a-Ce is still strongly correlated. 



a consequence of the behavior inuf (not shown): Within 
DMET(H-I), n/ is pinned at 1 for decreasing volume un- 
til V~17 A"^, at which point it increases and the system 
becomes mixed valent. 

Turning to the local magnetic moment, our approx- 
imations [neglect of spin orbit, intra-atomic exchange, 
and the 4/ crystal field splitting in Eq. (^)] have more 
serious implications for this quantity than others, and so 
we can provide only an estimate. Consistent with these 
approximations we take 



ifmcr'^ifni' a' 



) = 



ny/14 if mcr = m'a' 
d/91 if mcr ^ m'a' 



(20) 



such that the local magnetic moment becomes 



^) = Y.^{n^^r,A-n^f^if) = - (2/13) d, (21) 



indicating whether a local spin moment exists. Note that 
this quantity does not contain information about long- 
range magnetic order, aside from the fact that a finite 
moment would be required for such order. Also note that 
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(ml) is unlikely to vanish. Even if one just statistically 
distributes electrons with arbitrary i, m, and a, some 
sites will have electrons with the same spin and thus {mj,) 
will be finite, but it will be smaller than its maximal value 
obtained in the localized regime where d is minimal. 

The spin, orbital, and total angular momentum expec- 
tations can be expressed as (5'f/} — (3/4) (m^), (if/) = 
12(m^), and {Jff) = (51/4)(to^) due to the degeneracies 



in Eq. (20). Note that in the atomic limit (n/'^l, d^O) 



these expressions correctly give Sif = 1/2 and Lif =3, al- 
though (Jff) averages over the two spin-orbit multiplets. 
Our DMFT(QMC) result for (m^) at T = 0.054eV is also 
provided in Fig. |l^ (bottom dotted curve with open tri- 
angles), where this quantity is seen to drop by 5% from 
the 7 to the a volume. This may b^compared to 11% 
and 12% dropg for the DMFT(NCA)EI and impurity An- 
derson modelQ calculations, respectively, based on their 
values of 71/ and d = w{f'^). High-energy neutron scat- 
tering experiments observe single-ion magnetic response 
from 0.8 ±0.1 4/ electrons in the a phase, suggesting alaa 
that much of the local moment persists into that phase.EiJ 
Such high-energy or "fast" probes can detect a local mo- 
ment even if it appears screened out in "slower" measure- 
ments like magnetic susceptibility. Note that the, at first 
view unexpected, increase in the DMFT(QMC) (m^) for 
the smallest volumes in Fig. |lj only reflects this same 
behavior in n/ (Fig. p^ ). 

The persistence of a still robust (albeit slightly re- 
duced) local 4/ moment into the a phase as suggested, 
here supports the Kondo Volume Collapse scenario,L3 
in that the observed temperature-independent Pauli- 
likc paramagnetism of the a phase can then arise 
when the valence electrons .^'^pij'" out these local mo- 
ments. .0|i.'bi|tally polarizecOEErEJ and self-interaction 
correctedliZniflo LDA results suggest that the moment 
actually collapses to nothing in the a phase of Ce and its 
analog in Pr. However, these calculations really mea- 
sure spin and orbital polarization analogous to (wz), 
and therefore describe a loss of magnetic order in the 
a-like phases without providing information about the 
local moment itself. Indeed, there can be a local mo- 
ment (to^) Gven in the fully uncorrelated limit, as noted 
earlier, since (m^) = n/ — (2/13)dmax = "/(1^"//14) 
can be significant away from empty or filled bands. Note 
that one may have temperature-independent paramag- 
netism in the presence of local moments both if there is 
correlated Kondo screening of these moments, as noted 
above, as well as by Pauli's original one-electron process 
in which only electrons in states near the Fermi level are 
free to respond to the field. The latter must dominate as 
one approaches the uncorrelated V = limit. 



VI. SUMMARY AND DISCUSSION 

We have calculated thermodynamic, spectral, and 
other properties of Ce metal over a wide range of vol- 



umes and temperatures using the merger of the local 
density approximation and dynamical mean field theory 
(LDA+DMFT). The DMFT self energy was generated 
by rigorous quantum Monte Carlo (QMC) techniques, in- 
cluding a new, faster implementation that has facilitated 
lower-temperature results and is described in detail. Our 
LDA-fDMFT results provide a comprehensive picture of 
correlation effects in compressed Ce, and their fundamen- 
tal role in the first-order a-7 transition. First results of 
this effort have been published in Ref. ^ 

At large volume, we find a Hubbard split 4/ spec- 
trum, the associated local magnetic moment, and an en- 
tropy refiecting the degeneracy in the moment direction. 
This phase is well described by the Hubbard-I approx- 
imation and its energy but not its entropy also agrees 
with the polarized Hartree-Fock solution. As volume is 
reduced, a quasiparticle or Abrikosov-Suhl resonance be- 
gins to develop at the Fermi level in the vicinity of the 
a-7 transition, and the entropy starts to drop. At the 
same time, the 4/ double occupation grows whereas the 
number of 4/ electrons remains close to one. The tem- 
perature dependence of the quasiparticle peak is consis- 
tent with a significantly larger Kondo temperature in 
the a phase than in the 7 phase, and the parameter- 
free LDA-I-DMFT spectra are in good agreement with 
experiment for both a- and 7-Ce. In the range where the 
quasiparticle peak grows dramatically, the correlation en- 
ergy as a function of volume is seen to have a negative 
curvature. This leads to a growing shallowness in the 
total energy as temperature is reduced and is consistent 
with the first-order 01-7 transition within our error bars. 
Our results suggest that the temperature dependence of 
the transition pressure is primarily due to the entropy. 
Finally, if the volume is reduced below that of the ambi- 
ent a phase, the quasiparticle peak grows at the expense 
of the Hubbard side bands, yet these Hubbard side bands 
persist even at the smallest volumes considered. 

The Mott transitioril (MT) and Kondo volume 
collapseE3 (KVC) scenarios are based on the one-band 
Hubbard and the periodic (or more approximate im- 
purity) Anderson models as paradigms. The classi- 
fication of our results in terms of these two stan- 
dard theories requires distinguishing between the more 
general intecpretation of the MT in the many body. 



community]^ e.g., apphed to pach materials as V2O3, 
and the ideas of JohanssonEil and members of the 
local-density fu nC i tiopa l-jyHaim^ as applied to the /- 
electron metals .ll^llaO'EJMEjfiil In the former case, cor- 
related solutions of both model Hamiltonians show sim- 
ilar features at finite temperature such as persistence 
into the more weakly correlated regiipe of the local mo- 
ment and residual Hubbard splitting,E3 just as seen here 
for a-Ce. The similarity between the two models can 
be understood from the following consideration: The 
conduction-electrons in the periodic Anderson model are 
non-interacting. Thus, they only enter quadratically in 
the effective action and can be integrated out by a simple 
Gauss integration. This results in an effective one-band 
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model for the / electrons of the periodic Anderson model 
which can behavecj very much like the Hubbard model 
not only at finite temperature, but, depending on the 
choice of f-d hybridization, also at T = Q. 

One might try to distinguish between the two scenar- 
ios by whether the transition is caused by changes of 
the /-/ (MT) or /-valence (KVC) hybridization. But, 
since realistic calculations like the present include both, 
this distinction is rather problematic. Another difference 
can be addressed unambiguously, i.e., whether the low- 
temperature 7 phase has an Abrikosov-Suhl resonance 
(KVC) or not (MT). We observe the former, in agree- 
ment with experiment The energy scale of this 7-Ce 
Abrikosov Suhl resonance is very small such that we ob- 
tain a thermally smeared out structure instead of a sharp 
resonance. The smallness of the energy scale also implies 
that the effect on the total energy is very minor. Be- 
cause of this, the low-temperature energy (but not the en- 
tropy) of 7-Ce may also be adequately described by static 
mean field techniques like our HF calculation as well as 



a number of local-dens^ 
bitallii-|j 

lda,bNji 



rized LDA,BE 



uactional modifications: or- 
L3L-self-interaction corrected 
and LDA-f U.E^Eil These approximations 
have a frequency-independent (static) self-energy and 
provide a splitting of the 4/-band into two bands by an 
(artificial) symmetry breakins— While our HF calcula- 
tions as well as LDA-l-U workESi for Ce give a transition 
at too small volume, one may drive the onset of the sym- 
metry breaking closer to the volume of the a-7 transition 
by reducing the the 4/ Coulomb interaction Uf. Such a 
reduced interaction strength is naturally achieved within 
the orbitally polarized LDA calculations which omit U f 
and employ the weaker intra-atomic exchange interaction 
to achieve the symmetry breaking. 

A major point of debate between the KVC scenario and 
Johansson's interpretation of the MT picture is whether 
the a phase of Ce is strongly correlated (KVC) or not 
(MT). Our results suggest that a-Ce is strongly corre- 
lated with a three peak structure consisting of the two 
Hubbard peaks and central quasiparticle bands as in the 
KVC picture.^ In contraS|^t^&-MT madel as advocated 
by JohanssonEl and otherslJoO'BBEl predicts a single 
peak associated with uncorrelated, band-like / electrons. 
We do see a rapid increase in double occupation d across 
the transition, which is consistent with the delocalization 
ideas of this MT scenario. However, the actual value of 
d in the a phase is far from uncorrelated, although it 



indicates considerably less correlation than in the KVC 
pictured It appears that this perspective of the MT is 
motivated by the LDA results, and that if one were to 
fully take into account electronic correlations, one would 
also observe a correlated three-peak solution as in the 
Hubbard model. Ej This correlated solution would also 
have preformed local magnetic moments, which would 
be screened at small energies on the scale of the width 
of the Abrikosov-Suhl resonance, as in the KVC picture, 
whereas the uncorrelated MT solution does not develop 
such local moments. 

Since we find a strongly correlated a phase, the ques- 
tion remains why the structural and volume dependence 
of the total energy in the a-Ce regime is so extremely 
well described by normal pararpagHetic LDA and its gra- 
dient corrected improvements. PI Ml This point is-oiie of 
the st rO | iiges| . ysyiaenta advanced by JohanssorJij and 
othersBESyyOHBEl that a-Ce-like phases should be 
weakly correlated. A logical explanation would be that 
LDA may get the interactions between the quasiparticles 
correct but not their formation energy. The interactions 
are perhaps governed by the significant weight in the cen- 
tral Fermi-level peak which resembles the uncorrelated 
LDA solution, while the formation energy may involve 
the residual Hubbard sidebands in some way. Thus the 
still very significant correlation may provide only a con- 
stant contribution to the total energy in the a-Ce regime, 
so that the volume and structural dependence is still well 
represented. This would be consistent with the energy 
shift between a and 7 phases employed by Johansson et 
al.tZl in their LDA-based modeling of the transition. 
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